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Motivated by the various physical properties observed in #-(BEDT-TTF)2X, we study the 
ground state of extended Hubbard model on two-dimensional anisotropic triangular lattice 
at 1/4- filling with variational Monte Carlo method. It is shown that the nearest-neighbor 
Coulomb interaction enhances the charge fluctuation and it induces the anomalous state such as 
charge-ordered metallic state and the triplet next-nearest-neighbor /-wave superconductivity. 
We discuss the relation to the real materials and propose the unified view of the family of 
6»-(BEDT-TTF) 2 X. 
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Charge order (CO) is one of the important phenomena 
in strongly correlated electron systems as well as super- 
conductivity and magnetism. In a charge-ordered system, 
the cooperation or competition between charge and spin 
degrees of freedom cause the various interesting proper- 
ties as seen in manganese oxides, 1 ' 2 vanadium oxides, 3 
high-T c cuprates 4 and organic conductors. 5,6 Above all, 
the family of organic conductors 6>-(BEDT-TTF) 2 X (we 
denote #-(ET)2X in the following) is known to exhibit 
some anomalous COs 7 ' 8 and it is suggested that the su- 
perconductivity is induced by charge fluctuation in the 
vicinity of the CO state. 9 

6>-(ET) 2 X consists of alternating stack of X~ (closed 
shell) and ET 1 /^ (3/4- filled) layers and the latter con- 
tribute electrical conductivity. ET layer can be regarded 
as two-dimensional anisotropic triangular lattice with 
3/4-filled band. In a system away from half- filling, the 
charge degrees of freedom become important as well as 
the spin degrees of freedom. In addition, we have to take 
into account the long-range Coulomb interaction because 
the screening effect is rather weak in this system. To 
treat this problem, an extended Hubbard model which 
includes the nearest-neighbor Coulomb interaction V^ 
have been studied extensively. The mean-field study by 
Seo 10 successfully reproduces the horizontal stripe CO 
(insulator) in X=RbZn(SCN)4 confirmed by some exper- 
iments and supports the validity of the model. On the 
other hand, recent X-ray diffraction experiments have 
shown the existence of non-stripe type CO in some ma- 
terials. In X=RbZn(SCN)4, short-range CO with wave 
vector q = (1/3, k, 1 /4) is observed in a high temperature 
metallic region. 8 In X=CsCo(SCN)4, the coexistence of 
different COs, q 1 = (2/3, k, 1/3) and q 2 = (0,fc,l/2), 
is suggested. 7 These non-stripe type COs are not long- 
range order and only seen in a metallic region. It is con- 
sidered that the large value of nearest-neighbor Coulomb 
interaction Vij and the geometrical frustration cause the 
charge fluctuation and the anomolous COs. Theoreti- 
cally, several mean-field studies 11,12 for the extended 
Hubbard model suggest the existence of non-strype CO 
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Fig. 1. (a) Schematic representation of the arrangement of donor 
ET molocules. p and c denote vertical and diagonal directions, 
respectively, (b) Equivalent frustrated square lattice. 



with "3-fold" periodicity in the relevant parameter spaces 
for real materials. 

As for the superconductivity in X=T-3, its detailed 
property is still uncertain. Some theoretical studies dis- 
cuss the charge fluctuation mechanism caused by the 
nearest-neighbor Coulomb interaction V. In a square lat- 
tice, the large value of V enhances the charge fluctuation 
at (ir, 7r) and it leads to the checkerboard-type CO 13, 14 
and the d xy -wa,ve superconductivity. 15, 16 On the other 
hand, in a triangular lattice, the RPA study shows that 
V enhances the charge fluctuation at (2n/3, 2ir/3) and 
the triplet next-nearest-neighbor /-wave superconductiv- 
ity is induced. 17 This result is quite interesting since the 
triplet superconductivity is hard to be realized compared 
with the singlet one in general. However, the above stud- 
ies are limited to the mean-field level. Therefore, we need 
more accurate calculations which can treat the charge 
fluctuation more correctly. 

With these facts in mind, we study the following ex- 
tended Hubbard model in a two-dimensional anisotropic 
triangular lattice as the relevant model of 0-(ET)2X, 
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where denotes the summation over the nearest- 

neighbor sites. As shown in Fig. 1(a), the values of tij 
and Vij depend on their directions: t v and V v for the di- 
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Fig. 2. Three different CO patterns studied in this paper, (a) 
Vertical stripe and (b) Diagonal stripe state. Circles with up 
or down arrow denote charge-rich sites and open circles denote 
charge-poor sites, (c) The 3-fold state where black (gray) circles 
represent charge-rich A (charge-poor B and C) sublattices. 



agonal direction and t c and V c for the vertical direction, 
respectively. This lattice structure is topologically equiv- 
alent to a frustrated square lattice as shown in Fig. 1(b). 
We set U/tp = 10 and t c /t p — as the relevant values for 
X=CsCo(SCN) 4 (|*cApI < f) 18 and vaf y the valu es of 
V p and V c . To investigate the ground state of this Hamil- 
tonian, we use variational Monte Carlo (VMC) method. 
We introduce the following Jastrow-type trial wave func- 
tion, 



where 



I*) = PvPwPg |$) , 

P w = ^(p) Bi "i, 

P = pEw^i. 



(2) 

(3) 
(4) 
(5) 



Pg is a Gutzwiller projection operator which reduces the 
probability of double occupancy of electrons at the same 
site. P w and P v arc nearest-neighbor projection oper- 
ators which control the weight of the electron configu- 
ration at the nearest-neighbor sites. Summations over 
(p) and (c) correspond to the diagonal and vertical di- 
rections, respectively. The parameters g,w and v are 
variational parameters with < g,w,v < 1. If we set 
g = w = v = 1, is reduced to |3>) which is constructed 
from a mean-field solution of a certain Hamiltonian such 
as paramagnetic metal, CO and superconductivity. 

For the trial wave function of CO, we consider three 
different patterns as shown in Fig. 2. We call them (a) 
vertical stripe, (b) diagonal stripe and (c) 3-fold state, re- 
spectively. When V p 3> V c , the vertical stripe is expected 
because this pattern fully avoids the energy loss of Vp. 
For the same reason, the diagonal stripe is expected when 



Dashed circle represents the region where the superconductivity 
can appear. In Fig. 4, we show the energy comparison of three 
COs along the dashed arrow in the center of the phase diagram. 



Vp C V c . On the other hand, when V p ~ V c , the possibil- 
ity of non-stripe type CO exists since the stripe CO does 
not avoid both V p and V c simultaneously. Although we 
can consider other various types of non-stripe COs with 
long periodicity, we study the so-called "3-fold" state 
proposed by Mori 11 and discussed in a mean-field study 12 
as the first step. The 3-fold state consists of charge-rich 
A sublattice and charge-poor B and C sublattices, as 
shown in Fig. 2(c). This state is always metallic because 
the charge gap opens away from the Fermi energy. It is 
due to the fact that the 3-fold state is incommensurate 
for the triangular lattice with 1/4- filled band. 

For the trial wave function of superconductivity, we 
adopt the BCS-type wave function. The wave function 
with fixed electron numbers is used for calculational con- 
venience. The explicit form is given previously. 19 

Figure 3 shows the phase diagram in V p — V c plane 
obtained in the present VMC calculations. The vertical 
stripe is stable in the lower right-hand region (Vp 3> V c ) 
and the diagonal stripe is stable in the upper left-hand 
region (V p <C Vc), as expected. This result is consis- 
tent with the exact diagonalization study. 20 In the upper 
right-hand region, where the nearest-neighbor interac- 
tions are competing and large, the 3-fold metallic state is 
stable. This state is always metallic as mentioned before 
and have no spin order due to the geometrical frustra- 
tion. In the center of the phase diagram, on the other 
hand, the stripes and the 3-fold state have almost com- 
parable condensation energies as shown in Fig. 4. Inter- 
estingly, all of them are not insulating because they have 
semi-metallic Fermi surface with hole pockets and elec- 
tron pockets as shown in Fig. 5. These "charge-ordered 
metal" states are intermediate ones caused by the frus- 
tration of nearest-neighbor Coulomb interaction. 

Let us here compare our results with experiments. In 
X=CsCo(SCN)4, the coexistence of different short-range 
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Fig. 4. The relative energy compared to the 3-fold state. The 
horizontal axis is along the dashed arrow shown in Fig. 3. The 
system is in insulating state at both end of the axis and in semi 
metallic state at the center. 

CO patterns (the horizontal stripe and the 3-fold type) 
is observed in X-ray measurements. 7 This means that 
the stripes and the 3-fold state are energetically compet- 
ing with each other, as our result shows, and neither of 
them can become a long-range order at finite tempera- 
tures. For more accurate discussion, however, it will be 
necessary to consider the lattice degrees of freedom. The 
horizontal stripe observed in the experiments has higher 
energy than the diagonal one and does not appear in our 
phase diagram (Fig. 3). However, it is confirmed both ex- 
perimentally 8 and theoretically 10 ' 12 that the horizontal 
stripe is greatly stabilized in the #d-phase caused by the 
lattice distortion. Indeed, the horizontal stripe insulat- 
ing state is realized in X=RbZn(SCN)4 with structural 
phase transtion to ^-structure. 8 Thus we speculate that, 
in X=CsCo(SCN)4, the lattice distortion occurs locally 
and the horizontal stripe is stabilized, although it cannot 
be a long-range order since the 3-fold state has alomost 
the same energy and prevents the growth of horizontal 
stripe. 

Next, we discuss the possibility of superconductivity. 
As shown in Fig 3, we find that the novel type of super- 
conductivity has comparable condensation energy with 
those of the stripes and the 3-fold state in the lower left- 
hand region. Since the variational energies are so close 
to each other, it is rather difficult to determine which 
state is most stable. The pairing symmetry is the triplet 
next-nearest-neighbor (nnn) /-wave. The gap function of 
this symmetry, A(fe), changes its sign six times in a k- 
space and it corresponds to the nnn pairing in the real 
space. 19 The reason why the nnn /-wave is enhanced is as 
follows. When the nearest-neighbor Coulomb interaction 
is introduced, electrons avoid each other and nearest- 
neighbor pairing is suppressed. As a result, most stable 
pairing becomes possible on the nnn sites like the 3-fold 
state (Fig 2(c)). If the interaction is nearly isotropic, 
namely V p ~ V c , the 3-fold type charge fluctuation at 
q = (27r/3, 27r/3) is enhanced and nnn pairing is favored. 
The RPA study shows the same result. 17 This "charge- 
fluctuation-induced" superconductivity is a novel one in 




(a) (b) 

Fig. 5. Schematic representation of Fermi surfaces in the equiva- 
lent square lattice (Fig. 1(b)) for (a) the vertical stripe and (b) 
the 3-fold state, respectively. The black (gray) regions represent 
electron (hole) pockets. The dotted lines represent the original 
first Brillouin zone. 

contrast to the "spin-fluctuation-induced" d a ,2_ y 2-wave 
superconductivity in high-T c cuprates. 

In the previous paper, 19 we discussed the possibility 
of the nnn /-wave superconductivity in two-dimensional 
isotropic triangular lattice with electron density n — 2/3, 
having the superconductivity of Nao.3sCo02 • 1.3H 2 in 
mind. In that case, the 3-fold type CO is commensurate 
to the triangular lattice and greatly favored. Then the 
nnn /-wave has no chance to be stabilized since the con- 
densation energy of the 3-fold state is much larger. The 
1/iV expansion study by Motrunich and Lee also shows 
that the obtained T c is too small, although the most fa- 
vored pairing is the nnn /-wave. 21 In contrast, the sta- 
bility of the 3-fold state at n = 1/2 (1/4-filling) studied 
in this paper is much weaker than the case of n = 2/3 
because of the incommensurability. We can say that the 
stable 3-fold state is "melted" by the electron doping. In 
this condition, the nnn /-wave has comparable conden- 
sation energy and becomes the candidate of the ground 
state. It is possible that the superconductivity observed 
in X=l3 has nnn /-wave symmetry in the present mech- 
anism. Although we set t c /t p = in our calculation, the 
value of t c /t p in X=I 3 is estimated to be ~ 1.5. 22 t c /t p 
dependence of the stability of the nnn /-wave supercon- 
ductivity is an interesting future problem. 

Finally, we discuss an unified view of the family of 6- 
(ET) 2 X. When we replace the anion molecule X, the sym- 
metry of the system is kept unchanged and the change of 
the nearest-neighbor Coulomb interaction will be small. 
However, the hopping integrals change drastically and in 
particular the ratio t c /t p varies approximately from —0.5 
to 1.5. 18,22 We consider that the value of t c /t p plays im- 
portant roles and causes various physical properties. In 
the present paper, we study the case of t c /t p = and 
show the possibility of the coexistence of the stripes and 
the 3-fold state suggested in X=CsCo(SCN) 4 . For the 
case of t c /t p = —0.5, we obtain the result (not shown 
here) that the stripe COs are greatly favored and their 
region in the V p — V c diagram are enlarged compared with 
the case of t c /t p = 0. 23 Moreover, the region where the 
stripes and the 3-fold state are competing disappears. 
Combining with the mean-field study which includes the 
effect of lattice distortion, 12 we can explain the metal- 
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insulator transition with structural change observed in 
X=RbM'(SCN) 4 [M'=Co,Zn] and X=TlCo(SCN) 4 . For 
the case of t c /t p = 0.5, on the other hand, the stripes 
and the 3-fold state are suppressed and the paramagnetic 
metal state becomes dominant in the V p ~ V c region. 23 
Thus the value of t c /t p is crucial for the stability of CO 
and it will correspond to the horizontal axis of the exper- 
imental phase diagram of #-(ET) 2 X proposed by Mori et 
al. 5 

In summary, we have studied the extended Hubbard 
model on the two-dimensional anisotropic triangular lat- 
tice at 1/4- filling with variational Monte Carlo method. 
We have shown that for t c /t p = 0, the stripes and 
the 3-fold state are energetically competing with each 
other indicating the coexistence of different CO in 6- 
(ET) 2 CsCo(SCN) 4 . We also have shown that the 3- 
fold type charge fluctuation enhances the next-nearest- 
neighbor /-wave superconductivity and that it can be 
applicable to the superconductivity in #-(ET) 2 I 3 . The 
most important parameter is the value of t c /t p , which 
leads to a unified view of the family of 0-(ET) 2 X. 
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